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We present statistics of velocity fluctuations in both the Lagrangian and Eulerian frame for 
weakly driven two-dimensional turbulence. We find that simultaneous inverse energy and enstrophy 
ranges present in the Lagrangian and Eulerian Fourier spectra are not directly echoed in real-space 
moments of velocity difference. The spectral ranges, however, do line up very well with ratios of 
the real-space moments local exponents, indicating that though the real-space moments are not 
scaling "nicely" , the relative behavior of the velocity difference probability distribution functions is 
changing over very short ranges of length scales. Utilizing this technique we show that the ratios of 
the local exponents for Eulerian moments in weak two-dimensional turbulence behave in agreement 
with Kolmogorov predictions over the spectrally identified ranges. The Lagrangian local exponent 
ratios, however, behave in a different manner compared to their Eulerian counterparts, and deviate 
significantly from what would be expected from Kolmogorov predictions. 

PACS numbers: abc.123 



INTRODUCTION 

There are two reference frames that are normally con- 
sidered in turbulent fluids: the Eulerian frame and the 
Lagrangian frame [TJ The Eulerian frame of mo- 
tion is fixed to the laboratory frame where velocities, 
pressures and accelerations are fields fixed in space and 
varying with time {i.e., the velocity field u(x,i)). This 
frame has been used for many classical approaches to 
the turbulence problem, including the derivation of the 
Karman-Howarth equation and Kolmogorov scaling the- 
ory. The Lagrangian frame is fixed to fluid elements as 
they are advected by the turbulence. The position of the 
Lagrangian fluid element in time is denoted x(t), and 
the velocity and acceleration of the element are the flrst 
and second time derivative of the Lagrangian trajectory, 
respectively. The Lagrangian velocity at any time for 
any trajectory must equal its Eulerian counterpart {i.e., 
u{t) = u{x{t),t)). The Lagrangian frame is a useful one 
when considering the mixing of scalars [3j 131 [5], such 
as a dye, by turbulent motion. The difficulties in mak- 
ing experimental measurement of Lagrangian trajectories 
has made the characterization of Lagrangian turbulence 
much less common, with a few results reported for two- 
dimensional turbulence [6l [7] and for three-dimensional 
turbulent flows [SJ [S] . 

In either the Eulerian or Lagrangian frame, we would 
like to characterize the turbulent state by studying the 
statistical nature of the velocity flelds by measuring, for 
example, spectra and moments of velocity differences. 
We deflne the Eulerian ?i"^-order moment of longitudi- 
nal velocity difference as 

5(")(r)EE(((u(x + r)-u(x)).r)"), (1) 

and the Lagrangian n'^ order moment of velocity differ- 



ence as 

D(-Hr)^{\u{t + r)-u{t)n. (2) 

We would like to extract ranges of spatial scales r for 
the Eulerian frame and time scales r for the Lagrangian 
frame over which the velocity statistics exhibit scaling 
and measure the corresponding scaling exponents. These 
exponents are then compared with theoretical prediction, 
when such predictions exist. 

The moment characterization of the velocity fluctua- 
tions of turbulence is a standard procedure in fluid turbu- 
lence. Here we follow Frisch with respect to this stan- 
dard analysis paradigm. Velocity fluctuation moments 
have been experimentally obtained in three-dimensions 
for both the Eulerian and Lagrangian frames [H ITU] . 
Kolmogorov proved that for three-dimensional homoge- 
nous isotropic turbulence, S^^^{r) — — 4/5er, where e is 
the energy dissipation rate. Using dimensional analy- 
sis and the hypothesis of strict self-similarity, one can 
then show that S'*^"^(r) oc (er)"/'^. Subsequent measure- 
ments of fluctuations in three-dimensional turbulent flu- 
ids demonstrate that strict self-similarity does not hold 
in these systems, and that there are significant deviations 
from the expected n/3 scaling exponent, especially as r 
approaches the viscous dissipation scale and n becomes 
large |2] . This deviation is attributed to "intermittency" , 
a behavior in the spatial fluctuations of the velocity fields, 
characterized by bursts of activity. Intermittency is most 
commonly attributed to spatial fluctuations in the energy 
dissipation rate, e (note that e was assumed constant 
above) [2j . Accounting for this intermittency and the re- 
sultant adjustments to the Kolmogorov theory is still an 
active area of research. 

Intermittency is also a feature of velocity statistics in 
the Lagrangian frame of three-dimensional turbulence. 
Recent experimental [51 and numerical [IV measurements 
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of Lagrangian velocity difference statistics demonstrate 
stronger deviations from similarity than were displayed 
in the Eulerian statistics. Though there is no exact result 
for as there is for S'^"'\ a dimensional argument sug- 
gests that D'^^^t) cx (er)"/^ in the inertial energy range 
of scales. A more intermittent Lagrangian signal is per- 
haps not surprising; Kraichnan demonstrated |12j that 
even a non-intermittent Eulerian velocity field (albeit an 
unphysical one) could produce intermittent Lagrangian 
statistics. 

As in three-dimensional turbulence, there are exact re- 
lations and similarity results for Eulerian statistics in 
two-dimensional turbulence [T3j HH [151 US]- These re- 
sults apply to both the inertial energy range of two- 
dimensional turbulence and the inertial enstrophy range. 
Unlike three-dimensional turbulence, the Eulerian signal 
of the inertial energy range of two-dimensional turbu- 
lence is non-intermittent and self-similar, yielding scal- 
ing exponents for 5*^"-' (r) in agreement with the expected 
^n/3 [X21[T3]. Similarly, scaling exponents of ^^"^(r) for 
the Eulerian enstrophy range have the expected r" be- 
havior of the smooth velocity field at those scales. Al- 
though the scaling exponents in the enstrophy range are 
self-similar, the velocity difference moment may not be an 
adequate measurement tool to determine self-similarity 
for such smooth fields. This observation has resulted in 
the investigation of "inverse statistics" [T7]. Analysis of 
our data using the inverse-statistics approach will be pre- 
sented elsewhere. 

Whereas Eulerian results have been measured pre- 
viously for 2D turbulence |13l I15j . corresponding re- 
sults for Lagrangian velocity statistics have not been 
reported, perhaps owing to the difficulty in producing 
a two-dimensional system with significantly large scal- 
ing ranges for both energy and enstrophy. In three- 
dimensions, there is only a single inertial range in which 
the energy cascade range lies between the outer scale and 
the viscous scale. In two-dimensions there are two such 
ranges: an inverse-energy range analagous to the three di- 
mensional energy range, and the direct enstrophy range. 
One therefore needs inertial behavior and measurements 
over a larger range of scales in two-dimensional systems 
than is needed in three. Unfortunately, it is typical for 
two-dimensional systems to have small ranges, which re- 
sults in ambiguous scaling and a lack of self-consistency 
in how the ranges are defined. 

When the ranges in question are small and the mo- 
ments do not show unambigous scaling behavior, no 
strong conclusions about scaling exponents can be made. 
Weaker conclusions can be drawn, however, using tech- 
niques such as extended self-similarity (ESS) [H], which 
explores the shape of the probability distribution func- 
tions (PDFs) within a range of scales. If the PDFs ex- 
hibit a range over which their shape is either constant, or 
has a particular trend, it is possible to extrapolate to the 
case where one had a more extended range with a given 



scaling behavior. 

In this paper we report measurements of velocity dif- 
ference statistics in both the Eulerian and Lagrangian 
frames of two-dimensional turbulence. Moreover, we re- 
port these results for simultaneous energy and enstrophy 
ranges in our current system that is not significantly dif- 
ferent from prior experiments. Owing to the limited ex- 
tent of the ranges, we do not attempt to directly extract 
exponents of moments. Instead, we present ESS analy- 
sis of the local scaling exponents and demonstrate that 
the shape of the difference PDF's evolves very differently 
in the Eulerian and Lagrangian frames. If this shape 
evolution applies for larger ranges, we can conclude that 
the velocity difference statistics in the Eulerian frame for 
both the energy and enstrophy scaling ranges are self- 
similar, whereas in the Lagrangian frame they are not. 



EXPERIMENTAL 

Two-dimensional turbulence possessing simultaneous 
inverse energy and direct enstrophy cascades must be 
continuously forced [T3]. One method for such forcing, 
originally pioneered by Dolzhansky in 1979 [20], is to sub- 
ject a current carrying fluid to external magnetic fields. 
This method of forcing has since evolved into the strat- 
ified electromagnetic layer [21] which has become a com- 
mon systems for the study of 2D turbulence. 

The stratified electromagnetic layer consists of a dense 
salt water layer, typically 0.3 cm deep and 20 cm x 20 
cm square, underneath a less dense solution of roughly 
the same depth. A current is passed in the plane of the 
layer, and the layers are subject to a spatially varying 
magnetic field penetrating them vertically. The resultant 
Lorentz force drives fluid motion. The evolution of the 
stratified layer system is expected to approximate the 
forced/damped 2D Navier-Stokes equation. 
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supplemented by the incompressibility condition diUi — 
0, where u is the fluid velocity field, p is the density 
normalized pressure, F is the external electromagnetic 
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FIG. 1: Schematic illustration of the stratified-layer 2D tur- 
bulence experiment. 
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FIG. 2: Raw particle tracks obtained from the stratified layer 
over four consecutive frames (0.07 s). 



forcing, v is the fluid kinematic viscosity, and Einstein 
summation is used throughout. The hnear term with 
coefficient a represents the effects of frictional drag owing 
to the container bottom. Although such a description is 
broadly consistent with a number of experimental results, 
the detailed correspondence of the electromagnetic layer 
flow with the linearly-damped 2D Navier-Stokes equation 
has not been performed jT] [13] . 

In the set of experiments presented here, the system de- 
scribed above is modified slightly by replacing the lower 
layer of fluid with Fluorinert FC-75 and the upper layer 
by a dense salt water solution of 20% by mass NaCl with 
a small amount of liquid detergent added to lower sur- 
face tension and help with dissolution of tracer parti- 
cles. The Fluorinert has a density 1.8 times that of wa- 
ter with nearly the same viscosity, which allows for much 
stronger stratification than in the case of two salt-water 
solutions. It is also a strong dielectric so that only the 
upper salt-water layer is electromagnetically forced. Fi- 
nally, and perhaps most importantly, the Fluorinert and 
water are immiscible. These combined features allow the 
Fluorinert system to maintain stratification indefinitely, 
allowing the salt water layer to be driven harder than in 
previous systems PT] . 

The experimental apparatus, shown schematically in 
Fig.[l] consists of a 0.3 cm thick layer of salt-solution sus- 
pended over a 0.3 cm thick layer of Fluorinert. The layers 
are contained in a 40 cm x 20 cm box with reservoirs at 
each end across which copper electrodes are placed in 
the fluid. An alternating square-wave current with fre- 
quency 0.5 Hz and amplitude 0.75 amps is driven through 
the salt solution. A set of 1.27 cm diameter rare-earth 
magnets of approximately 0.7 T residual fleld strength 
are arranged with alternating field direction in a 20 cm 



X 20 cm square array with a period of 2.54 cm and ori- 
ented at 45° with respect to the current direction. The 
combination of the current and the magnetic field pro- 
duces a Kolmogorov-like forcing [TJ |5D] of alternating 
shear bands with the shear direction along y and peri- 
odicity Tinj « 1.8 cm in the x direction, which implies 
kinj = iT^jrinj =3.5 rad/cm. Using the layer depths of 
0.3 cm yields a « 0.125 Hz for the frictional coupling as- 
suming a simple linear shear in the Fluorinert. The salt 
solution upper layer has a viscosity around 1.15 that of 
water. 

To obtain Eulerian velocity fields and particle trajec- 
tory information, the upper salt-solution layer was seeded 
with polycrystalline powder with mean diameter 75 [im 
and density 0.98 gm/cc. Images of the particle fields, 
illuminated using several Xenon short-arc fiash lamps, 
were obtained with a 1280 x 1024 pixel CCD camera at 
a frame rate of 60 Hz. The velocity field was obtained 
from image pairs using particle tracking velocimetry de- 
rived from two earlier methods ^251 [55]. From a typical 
pair of images, of order 3 x 10* particle tracks were ob- 
tained and then interpolated to a 126 x 100 velocity field 
array. A typical field of raw particle tracks is shown in 
Fig.H 

From this raw data one can obtain particle trajectories 
in two ways: splice together raw particle tracks or use the 
interpolated Eulerian fields to generate particle trajecto- 
ries by solving the advection equation. The later tech- 
nique will be used here. The generated particle tracks 
solve the equation 

f = u(r(i),t), (4) 

where r(i) is the position of the tracer particle at time 
t, bicubic interpolation is used to approximate u at the 
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FIG. 3: Typical time-spliced real particle trajectory and a 
trajectory computed from dynamic velocity field data . 
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particle position, and fourth-order Runga Kutta is used 
to perform the time integration. Typical time-spliced 
particle trajectories and generated trajectories are shown 
in Fig. d 

The turbulent state that we obtain for the strati- 
fied layer system described above can be characterized 
by several of its Eulerian ensemble-averaged statistics. 
The average energy per unit mass E — is 8.4 

cm^/s^ (urms = 4.1 cm/s) and the average enstrophy 
il = w^„j^/2, where the vorticity w = V x u, is 51 s~^ 

{Wrnis = 10.1 S~^). 

RESULTS AND DISCUSSION 

Our system displays two different ranges corresponding 
to predominantly inverse-energy transfer [24, for scales 
larger than the injection scale and direct-enstrophy trans- 
fer for scales smaller than the injection scale. Because of 
the relatively limited range of scales available in the ex- 
periments, it is important to firmly establish the extent 
of the spatial and temporal scales for each transfer range. 
Empirically, this is best done utilizing spectra: we find 
that turbulent range transitions are sharper in Fourier 
space than in real space. The Eulerian and Lagrangian 
energy spectrum for data obtained from the electromag- 
netic cell are shown in Figs. |4^ and b, respectively. For 
the Eulerian spectrum, E{k) is calculated for individ- 
ual velocity fields using Ui{k)u*{k)/2 and averaged over 
velocity fields from the entire run. The Lagrangian spec- 
trum, E(f), is calculated from the Lagrangian velocity 
correlation {ui(t)ui(t+T)) using the Wiener-Kinchin the- 
orem. 

We begin by identifying the length scales over which 
there is an inverse-energy range. According to previous 
results [13l [HI [19], this range occurs at wave numbers 
smaller than the injection wave number ki„j = 27r/ri„j = 
3.5 rad/cm and is characterized by a spectral scaling of 
E{k) cx To identify this range, we plot in the inset 

the local spectral exponent, 

^ d{log{E{k)/E{k,^,))} 
d{logik/k„,,)} ■ 

The value —5/3 does not occur over any extended range. 
There is, however, a tight band of wave numbers smaller 
over which the values of x_e are close to —5/3, 
namely 1.65 cm~^ < k < 2.64 cm~^ (the dark bar in 
the inset plot). At wave numbers smaller than this range 
the spectral slope rises through zero and becomes pos- 
itive corresponding to the large-scale saturation peak, 
wave numbers just in excess of this range are associated 
with the injection scales. We take this tight band as 
the inverse-energy range; the choice of this tight band 
as a "range" will be justified in later analysis. A best 
fit power law to the range yields a spectral exponent of 
— 1.2. Taking the point at fc = 2.64 cm~^ and assuming 



E{k) = CQe^/^k~^/^ , where e is the scale-to-scale en- 
ergy transport rate (roughly the large-scale dissipation 
au^^^), yields a value for the Eulerian-Kolmogorov con- 
stant of Cq « 6, in agreement with earlier numerical and 
experimental work [T3J [TJ] . 

The same procedure is used to establish the direct- 
enstrophy range at wave numbers greater than kinj ■ The 
spectra in the direct-enstrophy range is expected to scale 
as E{k) cx k^^ in the absence of frictional dissipation [T^ . 
With frictional dissipation, as is the case in the electro- 
magnetic cell, the scaling exponent is expected to become 
more negative |15j . For wave numbers 4.5 cm^^ < fc < 13 
cm~^, indicated by a light bar in the inset, x has an 
average value of —5.7. This range is taken to be the 
direct-enstrophy range. Note that the scaling exponent 
is extremely steep when compared to the dissipation free 
expectation of —3, indicative of a large amount of exter- 
nal dissipation. The upper-limit of the direct-enstrophy 
range fc = 13 cm~^ corresponds to a length scale of 5 
mm, which is comparable to the fiuid layer depth of 3 
mm where one might expect 3D effects to become im- 
portant. 

Establishing the time scales of the inverse-energy and 
direct-enstrophy ranges in the Lagrangian frame is done 
in the same manner as for Eulerian statistics. The 
dimensional expectation for the energy range is that 
E{f) (X for frequencies smaller than the energy injec- 
tion frequency, fi^j = 1.66 Hz, which is about the eddy 
rotation frequency of an injection scale vortex. As in the 
Eulerian spectra, the Lagrangian spectra do not achieve 
a range over which the spectral scaling, measured by the 
exponent xl defined in the same way as xe and shown in 
the inset, achieves a constant —2 value. It also does not 
have as clearly defined a range as the Eulerian spectrum. 
The region marked by the dark bar corresponding to fre- 
quencies between 0.66 Hz and 1.28 Hz, characterized by 
a spectral slope increase from a —2 value to about — 1, 
should be the inverse-energy-range. Below this range the 
spectral exponent is less than 1 corresponding to large- 
scale saturation, and above is the abrupt spectral expo- 
nent change associated with /i„j . A best fit power law to 
this range yields an exponent of —1.6. Taking the point 
at / = 1.28 Hz and assuming E{f) = CQef~'^ yields a 
Lagrangian-Kolmogorov constant of Cq ~ 10. The en- 
strophy range, for which a simple dimensional prediction 
does not exist, is taken to be the range of frequencies 
spanning 1.87 Hz to 4.27 Hz, over which the spectral 
exponent is around —4.9. The frequency 4.27 Hz corre- 
sponds to a 0.23 s time scale, which is the same as the 
Lagrangian correlation time calculated by integration of 
the normalized two-point correlation. The steep spectral 
exponent is in agreement with earlier geophysical obser- 
vations |25j- 

We now consider Eulerian and Lagrangian velocity 
statistics as measured by n"^ order moments of velocity 
differences, i.e., structure functions, ^("^(r) and Z?(")(r), 
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FIG. 4; (a) Eulerian and (b) Lagrangian energy spectra. Insets show the behavior of the respective spectra's local scaling 
exponent. The proposed ranges of wave numbers (Eulerian) and frequencies (Lagrangian) over which there is an inverse-energy 
range is denoted by a dark bar in the insets. The proposed direct-enstrophy range is likewise indicated by a light bar. Best fit 
power laws over the suggested ranges are displayed by lines with the scaling exponents given. See the text for the exact ranges 
and power laws. 



respectively. Averages for the structure function calcula- 
tion are taken over ensembles of realizations, and spatial 
and temporal homogeneity of signals is assumed so that 
the moments do not depend on absolute position x or 
time t. For the inverse energy range of two-dimensional 
turbulence, we have dimensional predictions for scaling 
exponents: S^^'^r) cx (er)"/^ and i:>(")(r) cx (er)"/2. 
In the enstrophy range for Eulerian statistics, we ex- 
pect S'^")(r) cx r", indicating a smooth velocity field for 
which the linear order term in the Taylor expansion is 
dominant. For Lagrangian statistics, at small time scales 
where acceleration is approximately constant, we expect 
D^'^\t) (X r". Thus, we expect scaling in the enstrophy 
range to have exponents between n and n/2. 

We begin by looking at the second moments (n=2) 
which are displayed in Fig.|5] Inset in these plots are the 
local exponents for the second moments, defined as 

(„) ^ d{;og(5(")/<„J} 

- d{log{r/n^,)} ' 

(„) ^ d{;off(jj(")/<„j} 

- d{log{T/n^,)} • 

A scaling range where ^ has a constant value is not ob- 
served. This is unsurprising given the small extent of the 
ranges as identified by spectral analysis. 

The second moments are frequently used in lieu of 
spectral analysis to identify scaling ranges. This is done 
because, in the limit of long scaling ranges, the results 
should yield similar ranges (the Wiener-Kinchen theorem 
links the spectra with the second-moments through the 
correlation function ). For small ranges, however, simply 
assuming this correspondence is dangerous. If we deter- 



mine temporal and spatial scales utilizing these moments 
by, for example, bracketing ranges with the expected scal- 
ing exponents, we get inconsistent results. In particular, 
we expect the Eulerian second moment to scale as r^/"^ in 
the inverse energy range. Looking for the region where 

(2) 

assumes a value of 0.66 would yield a length scale of 
around 1.44 cm, which is significantly smaller than the 
injection scale of 1.8 cm, i.e., in the enstrophy range of 
length scales. A similar phenomenon happens for time 
scales of the Lagrangian second moment. The time scale 

(2) 

at which assumes the expected value of 1 is well 
within the spectrally identified enstrophy range. 

Given the lack of correspondence between real and 
spectral space results, it is interesting to ask if ranges 
identified spectrally are refiected in their real space coun- 
terparts. It is also valid to consider why we use the spec- 
tral method for identifying ranges in the first place since 
it was an empirical decision to do so. In the rest of this 
paper, we demonstrate that although the spectral ranges 
are not easily identified in the direct moments, they cor- 
respond directly to features of the PDF shape, that is, 
when we investigate velocity statistics using ESS. For the 
moment, simply assuming the spectrally identified ranges 
are correct (marked by a dark and light bar again), we 
get best fit power laws to the second moments that sig- 
nificantly deviate from expected behavior, see Fig. [Sj 

The Eulerian and Lagrangian velocity-difference mo- 
ments of order 1 < n < 6 are shown in Fig. [6] and the 
insets show the local exponents with dark and light bars 
to indicate the spectrally-identified scaling ranges. To fit 
exponents confidently to high-order moments, scaling is 
needed over an extensive range. Unsurprisingly, based on 
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FIG. 5: (a) Eulerian and (b) Lagrangian second-order structure functions. Insets show the behavior of the respective structure 
function local scaling exponents. The ranges of length (Eulerian) and time (Lagrangian) for the inverse-energy regime, deter- 
mined from spectra, are denoted by dark bars in the insets. The direct-enstrophy ranges, likewise determined from spectra, are 
indicated by light bars. Best fit power laws over the suggested ranges are displayed by lines with the scaling exponents given. 
See the text for the exact ranges and power laws. 




FIG. 6: (a) Eulerian and (b) Lagrangian n*'^-order structure functions normalized by u^^s- Insets show the behavior of the 
respective structure function local scaling exponents. The ranges of length (Eulerian) and time (Lagrangian) for the inverse 
energy regime, determined from spectra, are denoted by dark bars in the insets. The enstrophy-ranges, likewise determined 
from spectra, are indicated by light bars. 



the previous spectral analysis, no such ranges are visible 
at moments of any order. Although long ranges are not 
seen for either the Eulerian or Lagrangian data, the expo- 
nents of the moments behave in a very similar manner at 
different order, suggesting distinct changes in the behav- 
ior of the local exponents over the spectrally identified 
ranges. Because constant exponents are not observed in 
the measured range, we are forced to draw conclusions 
based on ESS analysis. 

We normalize the local exponents by the second-order 



local exponent and scale by the expected value of n/2, 
i.e., we use ESS. By doing this normalization, we are 
looking at the change in PDF shape over a range of scales. 
This is shown in Fig. [7] There is very different behav- 
ior for the Lagrangian and Eulerian analysis. For self- 
similarity to hold over a range (that is, an unchanging 
PDF shape), the value of the normalized local exponent 
must be unity, for all n, over that range, i.e., the higher 
moments are simply scaling as a power of the second mo- 
ment. Given these observations, it is possible to draw 
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FIG. 7: The structure functions local exponent for orders 1 < n < 6 normalized by n^{2)/2 for (a) the Eulerian frame and (b) 
the Lagrangian frame. The ranges of length (Eulerian) and time (Lagrangian) for the inverse energy regime, determined from 
spectra, are denoted by dark bars in the insets. The enstrophy-ranges, likewise determined from spectra, are indicated by light 
bars. A range with value of unity for all of the normalized local exponent indicates that the PDF shape is constant. 



some stronger conclusions. 

The spectrally-identified ranges, which seemed arbi- 
trary and did not clearly correspond to any particular 
behavior in the real-space statistics, become more trans- 
parent. The enstrophy range in the Eulerian statistics 
is characterized by a near unity grouping of all the mo- 
ments. At the upper-end of the Eulerian enstrophy range 
we see a sharp change in the behavior of the exponents 
where they begin to deviate from one. This deviation 
reaches a peak at the low-end of the inverse energy range 
and quickly collapses back to near unity. The conclu- 
sion we draw from this behavior is that, in our system, 
both the energy and enstrophy range Eulerian velocity 
difference statistics are tending to behave self-similarly 
(i.e., normalized exponents tending to unity). This is in 
agreement with earlier observations of the unnormalized 
quantities. [131 HH US] . We speculate, given previous nu- 
merical work with more extensive ranges, that the small 
deviation from unity that we observe, in the limit of a 
long range well removed from injection effects, should 
settle down to unity. The average value of the exponent 
ratio for both the Eulerian energy and enstrophy ranges 
for a range of n is displayed in Fig. [8] as open symbols. 
A line marked K41 (for Kolmogorov 1941) indicates self- 
similar scaling. 

The story is quite different for the Lagrangian results. 
The spectral ranges still find support, but in almost ex- 
actly the opposite manner: divergences from unity grow 
in the energy and enstrophy ranges while being reduced 
in the intermediate (injection) range (note that the ver- 
tical scales are different between Fig.[7|L and b). Indeed, 
for the energy range, the maximal deviation from unity 
occurs in the center of the range, collapsing back to unity 
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FIG. 8: The normalized local exponents for the velocity dif- 
ference moments averaged over the inverse-energy and direct 
enstrophy ranges in both the Eulerian and Lagrangian frame. 
Strict self-similarity over a given range would be indicated by 
collapse of the moments on the K41 (Kolmogrov 1941) line. 

only after the outer scale is reached. This allows us to 
speculate that, even in the limit of long ranges, the en- 
ergy range will continue to deviate from unity and dis- 
play a lack of similarity. Interpretation of what occurs 
in the enstrophy range is somewhat more complicated. 
There is a deviation from unity, but whether this is a 
residual result of a small range and the close proximity 
to the injection scale is uncertain. In the limit of a long 
range, the majority of the enstrophy range may approach 
unity. This would be unsurprising because the spectral 
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statistics indicate a very smooth enstrophy range in the 
Lagrangian frame. Therefore, we expect velocity differ- 
ence statistics to be dominated by first-order terms in the 
series expansion as in the Eulerian case. Without further 
measurements it is difficult to extrapolate based on this 
data. As with the Eulerian data, the average local ex- 
ponent over the energy and enstrophy ranges for a range 
of n is displayed in Fig. [8] as close symbols. Deviation 
from similarity is readily apparent when compared with 
the Eulerian statistics. 

The data at very large scales and times are not well 
enough converged to draw any conclusions about statis- 
tical behavior past the outer scale. The dissipative ranges 
(small r and r), however, are well converged and observ- 
ably different. The Eulerian range is characterized by an 
increasingly non-self-similar dissipation range, whereas 
the Lagrangian frame quickly collapses to a state of sim- 
ilarity. For the Eulerian frame, this is most likely in- 
dicative of a spatially inhomogenous viscous dissipation 
- a similar observation is found in three-dimensional tur- 
bulence [5]. For the Lagrangian frame, the collapse to 
similarity is surprising. It may be misleading, however. 
At these short times, linear order terms in a expansion 
are dominant, and therefore velocity difference may not 
be able to capture deviations from similarity well, much 
as is the case for the velocity difference statistics in the 
Eulerian enstrophy range. 



CONCLUSIONS 

We have considered velocity difference statistics in 
two-dimensional turbulence for both the Eulerian and 
Lagrangian frames. Eulerian statistics are investigated 
and compared with earlier results, and novel Lagrangian 
statistics are presented. As in prior experiments, the 
extent of the inverse-energy range and direct-enstrophy 
range are limited. The limited range has two major ef- 
fects: the scaling exponents are different in magnitude 
from the dimensional predictions, and there is no eas- 
ily identifiable correspondence between ranges in spectra 
and behavior of moments. These observations are not 
unexpected for systems without extended scaling ranges. 

Further analysis using ESS, however, yields a number 
of important insights. First, the results in the Eulerian 
frame are consistent with the broadly self-similar results 
found in prior experiments and simulations for velocty 
difference statistics, in spite of a lack of scaling. Second, 
the Lagrangian frame velocity difference statistics are not 
self-similar, a conclusion more strongly evident in the 
inverse-energy range than in the direct-enstrophy range. 
Finally, the dissipation ranges are not self-similar in the 
Eulerian frame, whereas they collapse to self-similar form 
for the Lagrangian frame. 

In addition to our main conclusions, we can provide 
two further insights from an analysis perspective. First is 



that careful identification of spectral ranges are echoed in 
the behavior of the moments local exponent ratios fairly 
precisely. This was somewhat surprising since it has been 
believed that moment ranges do not necessarily line up 
in detail with spectral ranges, in particular when those 
ranges are small. The data here indicate otherwise. Sec- 
ond, the local exponent ratio is effectively a derivative 
quantity. It measures the relative change in the behavior 
of the PDF's shape over a range of scales. Even when the 
moment ranges are small, the shape of the corresponding 
PDF changes rapidly. 
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